AFIT/GE/ENG/99M-18 


TERRAIN  BACKSCATTERING  COEFFICIENT 
GENERATOR 

THESIS 

Ricardo  Mediavilla,  Second  Lieutenant,  USAF 
AFIT/GE/ENG/99M- 1 8 


Approved  for  public  release;  distribution  unlimited 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Piihiir  rpnnrtinn  hiirrien  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
aatherin^and  mSninr  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  other  aspect  of  this 

??  lection  onnTorSn  indudi^^^^^  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  foi,«oformahon  Opera  ions  a^^^  121 5  Jefferson 

Davis  Highway.  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 


1 .  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE 

March  1999 _ 


4.  TITLE  AND  SUBTITLE 

TERRAIN  BACKSCATTERING  COEFFICIENT  GENERATOR 


6.  AUTHOR(S) 

Ricardo  Mediavilla,  Second  Lieutenant,  USAF 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
Air  Force  Institute  of  Technology 
2950  P  ST. 

WPAFB  OH  45433-7765 


3.  REPORT  TYPE  AND  DATES  COVERED 
Master’s  Thesis 


5.  FUNDING  NUMBERS 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFIT/GE/ENG/99M-18 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESSIEST 

Capt.  Kelce  S.  Wilson 
AFRL/SNAS 

2010  Fifth  Street  Bldg.  23,  Area  B 
Wright-Patterson  AFB,  OH  45433 


1 1 .  SUPPLEMENTARY  NOTES 

Advisor:  Dr.  Vittal  P.  Pyati 
Phone:785-3636  x4576 

Email :  V  ittal .  Py  ati@afit .  af .  mil _ _ 

12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 
Approved  for  public  release;  distribution  unlimited. 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words) 

Using  already  available  data,  a  backscattering  coefficient  generator  (BCG)  for  several  types  of  terrain  and  measurement 
conditions  is  developed.  The  types  of  terrain  are  1)  soils  and  rocks,  2)  trees,  3)  grasses,  4)  shrubs,  5)  short  vegetation,  6) 
roads,  7)  urban  areas,  8)  dry  snow,  and  9)  wet  snow.  These  data  sets  typically  cover  incidence  angles  ranging  between  0  and 
80  degrees  where  0  is  normal  to  the  terrain.  Measurement  conditions  are  defmed  by  1)  incidence  angle,  2)  wave 
polarizations  configuration  (HH,  HV,  or  VV),  and  3)  radar  band  (L,  C,  S,  X,  Ku,  Ka,  or  W). 

The  BCG  output  matches  very  closely  the  first  and  second  moments  of  the  published  data.  A  modified  chi-square  goodness  ol 
fit  test  at  a  0.05  significance  level  is  introduced  to  automatically  validate  the  generated  output  with  measured  data  histograms. 
The  BCG  reproduces  100%  of  measured  mean  backscattering  coefficient  (BC)  and  standard  deviation  values.  The  modified 
validation  test  failed  to  reject  the  BCG  output  as  representative  of  the  measured  data  for  72%  of  simulated  distributions 
having  a  database-source-point  number  (N)  greater  than  100.  Noted  discrepancies  can  be  attributed  to  sparse  BC  data 
histograms  and  BCG  limitations  at  incidence  angles  less  than  20  degrees  for  certain  terrain  t)q)es. 


14.  SUBJECT  TERMS 

Backscattering  Coefficient,  Ground  Clutter,  Chi-square, 


15.  NUMBER  OF  PAGES 
87 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRAC  i 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED _ UNCLASSIFIED  _ _ UL _ 

'  *  "  standard  Form  298  (Rev.  2-89)  (EG) 


UNCLASSIFIED 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. 


AFrr/GE/ENG/99M-18 


TERRAIN  BACKSCATTERING  COEEHCIENT  GENERATOR 


THESIS 


Presented  to  the  Faculty  of  the  Graduate  School  of  Engineering 
Of  the  Air  Force  Institute  of  Technology 
Air  University 

Air  Education  and  Training  Command 
In  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Master  of  Science  in  Electrical  Engineering 

Ricardo  Mediavilla,  BS 
Second  Lieutenant,  USAF 

March  1999 


Approved  for  public  release,  distribution  unlimited 


AFIT/GE/ENG/99M-18 


TERRAIN  BACKSCATTERING  COEFFICIENT  GENERATOR 


Ricardo  Mediavilla,  BS 
Second  Lieutenant,  USAF 


Approved: 


llJU.  F-  M _ 

Dr.  Vittal  P.  Pyati 
Committee  Chairman 

Maj  Michael  A.  Temple,  Ph.D. 

Committee  Member 


Jrew  J.^erzuoli 
^ittee  Member 


3  I'l 

Date 

Date 


; 


Date 


Preface 


This  research  proposes  a  simple  backscattering  coefficient  (o°)  generator  for 
several  radar  bands,  types  of  terrain,  and  wave  polarizations  for  angles  of  incidence  up  to 
80°.  The  generator  is  based  on  an  extensive  database  of  g°  measurements  [8].  This 
generator  was  first  proposed  and  successfully  used  by  Wilson  for  nine  terrain  types,  X- 
band  radar,  and  VV  polarization  [11]. 

This  research  expands  the  generator’s  capabilities  by  adding  six  radar  bands  and 
the  three  types  of  wave  polarizations  (HH,  HV  and  VV)  for  all  nine  terrain  types.  A 
method  to  automate  validation  of  the  generated  distributions  based  on  a  modified  chi- 
square  test  was  developed. 

Parallel  to  this  research  effort  is  the  work  of  O'Connor  (to  be  published).  The 
mentioned  research  successfully  applied  modified  sea  clutter  models  to  extend  the 
measured  a°  curves  to  the  grazing  angle  region  (0  >  80°).  The  "production"  G°  generator 
features  this  improved  grazing  angle  simulation  capability. 

I  want  to  thank  Capt.  Kelce  Wilson  for  the  invaluable  help  to  create  a  tool  that 
will  find  good  use  among  the  people  interested  in  backscattering  coefficient 
measurements. 


11 


Acknowledgments 


I  would  like  to  wholeheartedly  thank  my  thesis  committee  for  their  support,  ideas, 
and  suggestions  along  the  way:  Dr.  Vittal  Pyati,  Maj.  Mike  Temple,  Dr.  Andrew  Terzuoli, 
and  my  sponsor,  Capt.  Kelce  Wilson.  A  thanks  goes  to  all  supporting  personnel  at  the 
Air  Force  Institute  of  Technology.  These  individuals  made  possible  this  research. 
Thanks  to  all  friends  who  have  made  this  experience  a  good  one  every  day.  Special 
thanks  to  my  wife,  Paola,  for  her  comprehension  and  care  through  this  intense  journey. 
Blessed  be  our  Father  in  Heaven  for  all  his  help. 


Ricardo  Mediavilla 


iii 


Table  of  Contents 


Preface . 

Acknowledgements . 

List  of  Figures . 

Abstract . 

1.  Introduction . 1 

Background . 1 

Problem  Statement . 3 

Scope . 4 

Assumptions . 4 

General  Approach . 4 

Sequence  of  Presentation . 5 

n.  Literature  Review . 6 

Backscattering  Coefficient  (a°) . 6 

Scattering  Mechanisms . 8 

Dependence  of  a°  on  Terrain  Characteristics . 9 

Statistics  of  Ground  Clutter . 1 1 

Approximate  a°  Models . 13 

Summary . 15 

ITT  Backscattering  Coefficient  Generator  Development . 16 

Model  Formulation . 16 

Properties  of  0°  Generating  Function . 20 

Data  Quality  and  Pre-Processing . 21 

Validation  of  Results . 25 

IV.  Experimental  Results . 27 

Experimental  Data  Validation . 27 

Soils  and  Rocks . 30 

Trees . 31 

Grasses . 33 

Shrubs . 34 

Short  Vegetation . 35 

Roads . 36 


IV 


Dry  Snow . 36 

Wet  Snow . 37 

Summary  of  Results . 39 

V.  Conclusions  and  Recommendations . 40 

Appendix:  Simulated  Data  Sets . 42 

Bibliography . 73 

Vita . 76 


V 


List  of  Figures 


Figure 

1.  Data  Collection  Geometry . 8 

2.  Simplified  Facet  Scattering  Model . 8 

3.  Typical  Curves . 9 

4.  Array  of  Mean  a®  Values . 18 

5.  Array  of  s(0)  Values . 18 

6.  Total  Error  Array . "  20 

7.  a”  Generator  Output  Histogram . 21 

8.  Percentile  Plot  of  a°  Generator  Output  in  Figure  7 . 21 

9.  Partition  of  Data  Points  by  Even  and  Odd  Incidence  Angles . 22 

10.  Partition  of  Data  Points  as  in  Previous  Figure . 22 

1 1.  a"  Curve  for  Soils . 23 

12.  o°  and  s(0)  Curves  for  Trees,  X-Band,  HV  Polarization . 23 

13.  Spline  Interpolation  of  c°  Data  Points  for  Short  Vegetation . 24 

14.  Validation  Results  for  Short  Vegetation  (Ku-Band) . 27 

15.  Validation  Results  for  Grass  (S-Band) . 28 

16.  Chi-square  Test  Success  Rate  Distribution . 28 

17.  Normalized  Chi-square  Non-rejection  Rate . 29 

18(a).  Validation  Results  for  Short  Vegetation  (S-Band) . 30 

18(b).  Validation  Results  for  Grasses  (C-Band) . 30 

19(a-i).  Example  Validation  Results  for  Soils  and  Rocks . 3 1 


VI 


20(a-i).  Example  Validation  Results  for  Trees . 32 

21(a-i).  Example  Validation  Results  for  Grasses . 33 

22(a-i).  Example  Validation  Results  for  Shrubs . 34 

23(a-g).  Example  Validation  Results  for  Short  Vegetation . 35 

24(a-c).  Example  Validation  Results  for  Roads . 36 

25(a-i).  Example  Validation  Results  for  Dry  Snow . 36 

26(a-i).  Example  Validation  Results  for  Wet  Snow . 37 


vii 


AFIT/GE/ENG/99M-18 


Abstract 


Using  already  available  data,  a  backscattering  coefficient  (c°)  generator  for 
several  types  of  terrain  and  measurement  conditions  is  developed.  The  t5^es  of  terrain 
are  1)  soils  and  rocks,  2)  trees,  3)  grasses,  4)  shrubs,  5)  short  vegetation,  6)  roads,  7) 
urban  areas,  8)  dry  snow,  and  9)  wet  snow.  These  data  sets  typically  cover  incidence 
angles  ranging  between  0°  and  80“  where  0“  is  normal  to  the  terrain.  Measurement 
conditions  are  defined  by  1)  incidence  angle,  2)  wave  polarizations  configuration  (HH, 
HV,  or  VV),  and  3)  radar  band  (L,  C,  S,  X,  Ku,  Ka,  or  W).  The  o°  generator  output 
matches  very  closely  the  first  and  second  moments  of  the  published  data. 

A  modified  chi-square  goodness  of  fit  test  at  a  0.05  significance  level  is 
introduced  to  automatically  validate  the  generated  output  with  measured  data  histograms. 
Contrary  to  an  ordinary  chi-square  test  previously  applied,  the  modified  test  rejection  rate 
was  consistent  with  expected  behavior  of  simulated  and  measured  data  graphs  as  the 
number  of  source  data  points  (N)  increased.  The  modified  test  failed  to  reject  the 
proposed  generator  as  representative  of  the  measured  data  for  an  average  72%  of  all 
simulated  distributions  having  N  >  100.  Noted  discrepancies  can  be  attributed  to  sparse 
a°  data  histograms  and  model  limitations  in  the  0°  <  0  <  20°  for  certain  terrain  types. 

The  a°  generator  is  a  realistic  simulation  tool  of  terrain  backscattering  statistics 
since  it  portrays  terrain's  true  probabilistic  nature. 


TERRAIN  BACKSCATTERING  COEFFICIENT  GENERATOR 


I.  Introduction 


Back2round 

The  backscattering  coefficient  (a°)  is  a  measure  of  the  electromagnetic  energy 
scattered  from  a  distributed  target  in  the  receiver  direction.  In  general,  is  dependent 
upon  target  shape,  degree  of  roughness,  dielectric  constant,  viewing  geometry,  scattering 
center  location,  and  the  attributes  of  the  transmitted  and  scattered  electromagnetic  waves 
(wavelength,  polarization,  and  incidence  angle)  [8].  In  this  research,  the  backscattered 
energy  from  a  terrain  patch  of  definite  area  is  of  interest. 

The  scattering  centers  of  a  terrain  patch  are  randomly  distributed  and  vary  in 
position  with  frequency  and  incidence  angle  [8].  These  characteristics  give  G°  its  random 
quality.  The  reader  can  envision  such  an  effect  by  placing  some  glitter  on  a  piece  of 
paper.  When  illuminated,  it  will  become  evident  that  at  any  specific  "look-angle"  only 
certain  glitter  grains  flicker.  These  "scattering  centers"  on  the  piece  of  paper  contribute  to 
the  paper's  "a°".  These  a°  characteristics  justify  its  description  as  a  mean  value 
calculated  upon  several  measurements  at  an  incidence  angle,  frequency,  and  wave 
polarization.  These  will  be  referred  as  measurement  conditions.  The  value  of  G°  is  also 
influenced  by  seasonal  variations  of  vegetation  biomass  [8].  This  inherent  dependence  of 
G°  on  terrain  characteristics  has  found  generalized  use  as  indicator  of  crop  maturity, 
terrain  humidity,  and  other  uses. 
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The  potential  applications  of  o°  have  motivated  numerous  studies  and  papers 
dealing  with  its  measurement  and  theoretical  models.  Measurement  of  a°  has  been  a 
tough  task  because  of  the  randonmess  of  the  process  itself,  as  well,  as  the  calibration 
issues  involved.  Early  efforts  by  Ohio  State  University  and  other  research  groups  focused 
mainly  in  the  10°  <  0  <  80°  where  6  (incidence  angle)  is  measured  from  normal  incidence. 
Experiments  to  collect  c°  data  at  grazing  angles  (80°  <  0  <  90°)  are  very  scarce  [7]. 
Nevertheless,  the  measurement  of  a°  has  become  an  active  area  of  radar  research  because 
of  its  potential  military  and  environmental  applications.  Military  applications  arising 
from  the  theory  and  measurements  of  a°  focus  in  the  development  of  accurate  target 
recognition  techniques.  Environment  sensing  applications  are  being  applied  to  optimize 
the  use  of  natural  resources. 

In  this  research,  the  measured  o°  data  source  is  the  Handbook  of  Radar  Scattering 
Statistics  for  Terrain  [8].  This  book  features  a  comprehensive  database  of  a° 
measurements  for  several  t5q)es  of  terrain  and  measurement  conditions.  These  data 
measurements  were  collected  from  several  studies  spanning  several  yeajrs  and  intents. 
The  database  portrays  the  statistical  properties  of  ct°  measurements  as  a  mean  a°  and 
associated  s(0).  Frequencies  of  occurrence  histograms  describe  the  distribution  of  the  a° 
values.  The  classifications  for  terrain  are:  (1)  soil  and  rock  surfaces,  (2)  trees,  (3)  grasses, 
(4)  shrubs,  (5)  short  vegetation,  (6)  roads,  (7)  urban  areas,  (8)  dry  snow,  and  (9)  wet 
snow.  The  conditions  for  data  inclusion  into  the  database  are  comprehensive  and  well 
explained  in  Chapter  V  of  the  mentioned  reference  [8:79]. 
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This  research  proposes  a  a°  generator  that  is  constructed  upon  the  compiled  data 
in  [8].  The  user  of  the  (regenerator  must  specify  (1)  type  of  terrain,  (2)  wave  polarization, 
(3)  incidence  angle,  (4)  radar  band,  and  (5)  a  f/(0,l)  random  variable.  These  will  be 
called  generator  input  set  (GIS).  The  generator  produces  a  value  of  <j°  in  decibels.  Note 
that  to  obtain  the  published  mean  a°  and  s(e)  for  a  GIS,  sufficient  a°  values  must  be 
generated. 

In  his  doctoral  dissertation,  Wilson  first  proposed  and  used  a  prototype  a° 
generator  for  all  types  of  terrain  listed  above  but  limited  to  X-band  and  VV  polarization 
[11].  This  research  applies  previously  established  concepts  to  expand  the  generator's 
capabilities  to  include  nine  terrain  types  (described  above),  six  additional  radar  bands, 
and  the  full  set  of  linear  wave  polarizations,  namely  horizontal-horizontal  (HH), 
horizontal-vertical  (HV),  and  vertical-vertical  (VV).  The  generator  output  is  tested  for 
goodness-of-fit  against  the  measured  data  available  using  a  chi-square  hypothesis  test 
with  encouraging  results. 

Problem  Statement 

Develop  a  backscattering  coefficient  (a°)  generator  that  will  emulate  measured  a° 
data.  The  generator  should  be  designed  with  a  minimum  1°  angular  resolution.  The 
generator  output  statistics  should  match  the  first  and  second  moments  of  the  published 
data  in  [8]  when  a  sufficient  number  of  samples  are  generated.  The  generator  should 
output  a°  values  for  all  three  pairs  of  linear  polarizations,  seven  radar  bands,  and  nine 
types  of  terrain.  Generated  ouput  distributions  should  be  similar  to  the  frequency  of 
occurrence  histograms  of  measured  data.  Validate  generated  distributions  using  a  chi- 
square  goodness-of-fit  test  for  all  incidence  angles  where  data  exists.  Generator  should 
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be  autonomous  such  that  new  a°  data  and  generated  histogram  validation  may  be  done 
with  minimum  user  interaction. 

Scope 

Limit  data  to  the  aforementioned  nine  terrain  types,  seven  radar  bands,  and  the  full 
set  of  linear  polarizations,  as  tabulated  in  [8].  The  a°  generator  algorithm  should  be 
coded  in  FORTRAN. 

Assumptions 

It  is  assumed  that  the  database  in  [7]  is  accurate. 

General  Approach 

The  c°  generating  function  is  defined  by 

a  =p{-ln{uT)  +  r  (1) 

where  P  is  an  amplitude  term,  a  is  a  skewness  parameter,  7  is  a  distribution  adjustment 
quantity,  and  U  is  a  uniform  17(0,1)  random  variable.  The  mean  of  the  function  over  an 
ensemble  of  f/(0,l)  random  numbers  in  (1)  represents  the  mean  a°,  while  its  standard 
deviation  represents  the  standard  deviation  of  the  mean  a°,  namely  s(6). 

The  mean  <f  and  s(0)  values  contained  in  [8:119-356]  are  typically  provided  in 
angular  increments  of  5  and  10  degrees.  A  method  of  interpolation  by  cubic  splines  is 
applied  to  the  data  to  achieve  the  required  1°  angular  resolution.  Cubic  spline 
interpolation  was  chosen  on  the  basis  of  smoothness  and  its  inherent  ability  to  connect  all 
measured  data  points. 

The  a,  P,  and  /values  are  selected  such  that  all  interpolated  values  of  the  mean  <f 
and  s(0)  sequences  can  be  generated  by  (1).  An  optimization  algorithm  is  used  to 
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determine  which  a-fi  combination  (for  a  fixed  /  value)  generates  the  desired  a°  and  s(0) 
with  minimum  error.  The  generated  distributions  are  tested  for  goodness-of-fit  using  a 
chi-square  test.  If  a  generated  distribution  is  rejected,  parameter  yis  either  increased  or 
decreased  to  adjust  the  output  histogram  of  (1)  such  that  it  fits  the  measured  data 
histograms  satisfactorily.  Manipulation  of  y  has  some  practical  limitations  that  are 
discussed  in  Chapter  IV. 

Sequence  of  Presentation 

A  literature  review  about  and  its  statistical  models  and  a  brief  description  of 
some  of  its  applications  take  place  in  Chapter  n.  A  description  of  Equation  (1)  will  be 
given  in  Chapter  HI.  The  algorithms  used  to  find  the  appropriate  a-j3  sets  are  developed 
and  explained  in  this  chapter.  Presentation  and  analysis  of  the  results  follow  in  Chapter 
IV.  To  end  this  journey,  conclusions  and  recommendations  are  given  in  Chapter  V. 
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II.  Literature  Review 


In  the  early  years  of  radar  development,  it  became  evident  that  unwanted  clutter 
originating  from  electromagnetic  waves  echoing  back  from  terrain  and  sea  needed  to  be 
accounted  for  in  radar  system  design.  Efforts  to  quantify  clutter  and  its  characteristics 
support  radar  systems  engineering  and  more  recently,  aid  in  the  monitoring  of  terrestrial 
phenomena. 

Backscattering  Coefficient  (a”) 

The  radar  range  equation  describes  the  relationship  between  the  transmitted  power 
Pt  to  a  point  target  at  a  range  R  and  the  power  Pr  detected  by  the  radar  receiver  [7] 


P.G, 

47cR^ 


(2) 


The  quantities  G,  and  Ae  represent  the  transmitting  antenna  gain  and  the  receiving 
antenna’s  effective  aperture  area,  respectively.  The  quantity  Orcs  is  the  target's  radar  cross 
section  (RCS)  in  units  of  meter  squared  (m^). 

The  RCS  of  a  target  depends  on  radar  system  characteristics  such  as  direction  of 
illumination,  wavelength,  and  polarization  [7].  Also,  a  target's  RCS  depend  on  the 
complex  permmitivity,  degree  of  roughness,  and  viewing  geometry.  In  mathematical 
terms,  the  RCS  is  defined  as  [8] 


r 


’  A 


V 


(3) 


where  the  subscript  p  and  q  refer  to  the  polarization  of  the  scattered  and  incident  waves, 
respectively.  The  quantity  R  is  the  radar  to  target  distance.  and  S'^  are  the  scattered 


and  incident  power  densities,  respectively.  The  scattered  power  density  (numerator  of 
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equation)  varies  as  rendering  Equation  (3)  independent  of  range  [7].  The  limit 
notation  denotes  that  the  RCS  is  a  far  field  quantity.  The  set  of  linear  polarizations  (HH, 
HV,  and  VV)  is  used  to  be  consistent  with  the  convention  in  [8]. 

By  using  the  relationship 

ri 


where  t]  is  the  air  impedance  and  E  is  the  electric  field  intensity,  we  have 


f 

K 

2> 

4n-R^- 

p 

Tpi 

r 

[  \ 

y 

(4b) 


where  p  and  q  are  substituted  by  HH,  HV,  or  VV  as  appropriate.  Note  the  reciprocity  in 
the  terrain  backscattering  process  for  HV  and  VH  polarizations  [8] 

<  =  <H  (5) 

A  more  convenient  measure  of  RCS  for  distributed  targets  is  the  backscattering 
coefficient,  cf,  defined  as  [8] 


o  ^RCS 

po-  ^ 


(6) 


where  A  is  the  area  covered  by  the  antenna's  footprint  on  a  reasonably  flat  surface  [10]. 
This  causes  a°  to  be  independent  of  the  radar's  geometric  parameters  (pulse  width, 
beamwidth,  etc.)  [7].  In  this  research,  the  definition  of  c°  as  in  Equation  (6)  will  be  used 
(Figure  1). 
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Figure  1.  Data  collection  geometry. 


Scattering  Mechanisms 

Several  theories  attempt  to  explain  the  interaction  of  electromagnetic  waves  and 
terrain.  Although  a  thorough  discussion  of  these  is  out  of  the  scope  of  this  document,  the 
facet  theory  of  scattering  will  be  briefly  described. 

This  theory  models  an  object  by  a 
number  of  facets  tangential  to  the  object's 
surface.  The  centers  of  the  facets  constitute 
scattering  centers.  The  total  scattering  is 
computed  by  determining  the  re-radiation 
patterns  and  slope  of  each  facet.  Figure  2 
Rgure  2.  Simplified  facet  scattering  model.  displays  a  highly  simplified  facet  model  for  a 

tree  canopy.  The  facets  can  be  large  or  small  relative  to  wavelength.  As  the  wavelength 
increases,  a  given  facet  will  change  from  large  to  small.  The  overall  effect  of  the  increase 
in  wavelength  is  the  apparent  smoothening  of  the  surface.  When  the  transmitted  signal  is 
narrow-band  and  coherent,  the  received  signal  will  be  the  phasor  sum  of  each  facet  (or 
scatterer)  contribution  [7]. 
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Dependence  of  a”  on  Terrain  Characteristics 


The  a°  values  for  a  distributed  target  are  determined  by  its  electromagnetic  and 
physical  attributes.  The  of  a  co-polarized  scatter  (HH  and  VV)  exhibits  three  distinct 
regions,  namely  "quasy-specular",  "plateau"  region,  and  a  "shadow"  region  (Figure  3).  In 

general,  the  (f  of  a 
cross-polarized  scatter  (HV)  does 
not  exhibit  the  same  distinct 
regions  as  the  co-polarized  scatter 
[7].  The  most  important  details 
about  cr°  response  based  on  terrain 
characteristics  will  now  be 
summarized.  This  summary  will 
follow  [8]  closely. 

For  smooth  surfaces,  the  exhibits  a  strong  return  in  the  quasi-specular  region 
that  decays  rapidly  as  the  angle  of  incidence,  0,  increases.  For  rough  surfaces,  the  change 
in  a°  as  6  increases  is  of  much  smaller  magnitude.  Generally,  the  ct°  curve  exhibits  a 
slight  negative  slope  [7]. 

Terrain  moisture  greatly  influences  the  amount  of  scattered  energy.  For  a 
distributed  target,  the  value  of  increases  as  the  moisture  content  increases  and 
viceversa  [8].  This  is  a  consequence  of  the  larger  microwave  dielectric  constant  of  water 
relative  to  that  of  the  holding  surface.  This  property  has  been  useful  as  a  terrain  humidity 
indicator.  For  instance,  a  study  at  the  Freiburg  test  site,  Upper  Rhine  Valley,  Germany 


Figure  3.  Typical Curves  [8]. 
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showed  strong  positive  correlation  between  o°  and  soil  moisture  content  at  microwave 
frequencies  [6].  A  satellite-based  system  was  used  to  measure  c°. 

For  the  most  part,  surfaces  covered  by  vegetation  have  a  (f  response  mainly  due 
to  the  canopy,  and  to  a  lesser  degree,  to  the  underlying  soil.  As  the  vegetation  canopy 
thins,  the  a°  contribution  from  the  soil  increases,  changing  a°  [7].  The  thinning  of 
vegetation  decreases  the  a°  response  of  a  distributed  target.  A  case  in  point  is  a  study  of 
Les  Landes  Forest,  France  (L  band)  [2].  This  study  showed  that  forest  thinning 
contributed  to  cf  reduction  of  up  to  4  dB,  most  noticeably  for  HV  polarization.  Thinning 
involves  decreasing  tree  population,  which  induces  a  change  in  mean  tree  height  and  a 
change  in  the  trunk  diameter  distribution.  Also,  thinning  causes  canopy  structure 
modifications  such  as  rearrangement  of  tree's  crowns  and  the  clearing  of  vegetation  above 
ground  level.  The  authors  theorize  that  this  could  be  responsible  for  the  g°  decrease. 

Snow  exhibits  a  c°  signature  that  changes  with  the  observation  wavelength.  If  the 
wavelength  is  large  compared  to  the  internal  structure  of  the  snow  layer  (such  as  L  band), 
the  a°  will  be  largely  from  the  underlying  soil.  This  changes  as  the  wavelength  becomes 
smaller  (Ka  band).  The  a°  response  will  resemble  that  of  a  rough  surface  since  the  snow 
structure  becomes  the  main  source  of  backscatter.  Moisture  also  has  a  significant  impact 
on  snow's  a°.  For  two  samples  of  snow  of  differing  moisture  content,  the  a°  response 
will  be  stronger  for  the  drier  sample  as  the  wavelength  increases.  This  difference  is 
greatest  at  small  wavelengths  (Ka  band)  [8]. 

Urban  areas  are  hard  to  describe  using  g°  since  many  structures  in  such  an 
environment  have  deterministic  shapes  that  can  behave  as  strong  point  targets  [7].  These 
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point  targets  scatter  energy  that  dominates  over  the  total  scattered  energy  from  all  other 
scatterers  present.  The  theory  to  handle  this  situation  is  different  to  that  of  other  terrain 
types.  For  example,  a  recent  study  of  the  urban  characteristics  of  Sydney,  Australia 
computed  theoretical  a°  signatures  from  urban  areas  using  deterministic  methods  [4]. 
This  investigation  identified  single  bounce  scattering  from  surfaces  such  as  roofs,  double 
bounce  from  ground-wall  structures,  and  other  metallic  elements  as  principal  radar 
scattering  contributors.  A  better  understanding  of  urban  area  scattering  will  enable 
monitoring  and  precise  classification  of  such  areas. 

Other  method  to  quantify  scattering  from  urban  areas  is  to  measure  the  c° 
response  of  an  urban  conglomerate  and  average  the  return  over  wide  areas.  Then,  the 
scanned  areas  could  be  classified  as  residential,  commercial,  or  industrial  [7]. 

Statistics  of  Ground  Clutter 

Most  clutter  theories  assume  that  a  patch  of  terrain  is  large  enough  to  contain  a 
significant  number  of  randomly  localized  scatterers  [7].  These  scatterers  produce  a  signal 
at  the  receiver  whose  envelope  behaves  like  random  noise  due  to  fading  effects.  (Fading 
refers  to  the  wide  variations  in  the  received  signal’s  envelope  due  to  phase-shift  changes; 
in  turn,  these  are  caused  by  look-angle  changes  relative  to  the  illuminated  area  due  to 
platform  movement.)  Therefore,  the  statistics  of  random  noise  can  be  used  to  model  the 
scattered  signal  being  received,  namely,  the  Rayleigh  clutter  model  [7]. 
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The  application  of  the  Rayleigh  clutter  model  is  based  on  a  set  of  assumptions  discussed 
in  [8]  and  reviewed  here  for  convenience: 

(1)  It  is  assumed  that  the  individual  scatterers  are  independent  of  each 
other.  Any  interactions  between  adjacent  scatterers  is  ignored.  Hence,  the  total 
instantaneous  electric  field  is  expressed  as  the  sum  of  each  scattered  electric  field. 

(2)  The  maximum  extent  of  the  target  under  scrutiny  is  much  smaller  than 
the  mean  range  from  the  antenna  to  the  target.  The  antenna  is  assumed  to  have 
uniform  gain  across  the  target  area. 

(3)  As  stated  before,  a  distributed  target  contains  a  large  number  of 
scatterers.  The  amplitude  and  phases  of  the  scattered  electric  fields  are  considered 
independent.  This  last  condition  is  satisfied  by  having  the  scatterers  randomly 
distributed  across  the  target's  area.  At  least  ten  scatterers  are  necessary  to  ensure 
that  the  orthogonal  components  of  the  electric  field  are  normally  distributed  [8]. 

(4)  The  phase  must  be  randomly  distributed  between  [0,  2k].  This  can  be 
accomplished  by  choosing  the  target's  extent  to  be  several  wavelengths  across. 
The  target  is  also  assumed  to  have  randomly  located  scatterers. 

(5)  No  single  scatterering  center  produces  an  electric  field  intensity  that 
dominates  the  total  scattered  field  from  the  distributed  target. 
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These  assumptions  lead  to  the  Rayleigh  clutter  model: 


/W  =  ^  OS0<27t 


where  is  the  probability  density  function  (PDF)  of  the  total  electric  field  magnitude 
before  detection  at  the  receiver.  The  second  equation, /(0)  is  the  PDF  for  the  phase,  while 
Ee  denotes  the  average  value  of  Eg.  The  variance  of  Ee  is  denoted  by  s(6)  and  the 

standard  deviation  ,  s(0),  is  just  the  square  root  of  the  variance. 

The  Rayleigh  clutter  model  may  not  be  applicable  to  forest  or  urban  scenes  where 
some  scatterers  may  dominate  the  return  at  certain  frequencies  and  incidence  angles  [8]. 
The  Ricean  PDF  is  used  when  one  or  a  few  strong  scatterers  exist  within  the  observed 
target.  Other  PDF  models  that  are  often  used  include  the  Weibull  and  log-normal. 
Approximate  a°  Models 

In  this  section,  the  various  theoretical  and  empirical  models  published  by  several 
a°  studies  will  be  explained. 

The  following  set  of  equations  describes  behavior  of  a°  of  the  combined 
observations  of  space  station  Skylab's  radiometer-scatterometer  (RADSCAT)  and 
University  of  Kansas'  truck  mounted  microwave  active  spectrometer  (MAS). 
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The  model  is  representative  of  summer  conditions  averaged  all  over  North  America: 

(j"Xf,0)=A  +  B-9  +  C-f  +  D-f-e  2O"<0<7O° 

cr:,{f,0)  =  M{0)  ^N{0)-f  e  =  0°,  10° 

where /is  the  frequency  and  0  is  the  incidence  angle.  The  indicated  range  of  incidence 
angles  for  both  equations  determines  the  values  of  the  A,  B,  C,  D,  M,  and  N  coefficients. 
The  loading  parameters  for  both  equations  can  be  found  in  [8].  The  a°  value  for  an 
incidence  angle  0°  <  0  <  10“  and  10°  <  0  <  20°  is  found  by  interpolating  the  resulting 
curve. 

Several  studies  have  investigated  the  potential  applications  of  radar  for  moisture 
sensing  of  vegetation  covered  soils.  The  following  model  represents  the  a°  response  of  a 
soil  surface  covered  by  vegetation  [9]: 

a  =T  -cr^  +o-^, 

where  is  the  two-way  attenuation  effect  due  to  the  vegetation  layer,  <f^  is  the 
backscattering  contribution  of  the  soil  underneath,  cj°<,v  is  the  direct  backscatter 
contribution  of  the  vegetation  cover,  and  o°jn,  represents  multiple  scattering  due  to  soil- 
vegetation  canopy  interactions. 

The  last  model  for  (t°  to  be  reviewed  models  the  ct°  measured  over  bare  soils, 
taking  into  account  soil  moisture  and  incidence  angle  [1]: 

cr"  =  C|  +  Cj  cos{0j^  +  dm^ 

where  C,,  C,,  and  Cj  correspond  to  the  backscattering  coefficient  of  dry  soil  for  0  =  90°, 
the  total  dynamics  of  soil  response  over  the  0°  <  0  <  60°  range,  and  a  shape  factor, 
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in  volumetric  soil  moisture.  Other  important  details  of  the  formulation  are  found  in  the 
reference  article. 

Summary 

The  papers  mentioned  in  this  chapter  represent  a  small  sample  of  the  available 
literature  related  to  (f  measurements.  Several  efforts  are  directed  towards  remote  sensing 
applications  that  seek  to  optimize  the  use  of  the  limited  earth  resources.  Also,  military 
applications  arising  from  the  theory  and  measurements  of  a”  logically  follows.  Target 
recognition  technologies  seek  to  precisely  identify  targets  camouflaged  on  the  ground. 
Hence,  an  understanding  of  ground  clutter,  as  quantified  by  (7°,  and  of  the  scattering 
mechanisms  is  crucial  to  create  effective  target  acquisition  algorithms.  In  the  next 
chapter,  a  method  of  simulating  <f  values  based  on  measured  data  is  introduced.  In 
contrast  to  some  models  described  in  this  chapter,  the  G°  generator  behaves  in  a  statistical 
manner,  yielding  histograms  similar  to  those  of  measured  data. 
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III.  Backscattering  Coefficient  Generator  Development 


The  developmental  aspects  of  the  a°  generator  are  explained  in  this  chapter.  The 
main  idea  behind  the  a°  generator  is  to  determine  a  functional  form  that  can  simulate  the 
first  and  second  moments  of  the  backscattering  process  for  various  terrain  types  in  a 
statistical  way.  By  "statistical"  it  is  meant  that  the  generator  returns  a  a“  number  that  is 
different  every  time  it  produces  such  a  value.  Further,  the  frequency  of  occurrence 
histograms  of  the  generated  a°  values  must  be  similar  to  those  of  measured  histograms 
when  the  generator  produces  enough  samples. 

Model  Formulation 

The  following  equation  is  the  functional  form  described  in  (1)  and  repeated  here 
for  convenience: 

cr  =fi-(-\n{U)")  +  r  (8) 

The  first  step  in  this  process  is  to  generate  an  array  (or  surface)  of  values  that  for 
possible  combinations  yields  the  required  mean  a"  and  s(9)  values.  Each  surface 

is  constructed  for  a  fixed  y  value.  To  illustrate,  a  typical  set  of  values  follows: 


a  =0.06,  0.061,0.062, ...  ,2.0 

(8a) 

P=  0.01, 0.011,  0.012,  0.013, ... ,  14.0 

(8b) 

y=0.01 

(8c) 

For  each  combination  of  a,  p,  and  y,  an  average  is  calculated  upon  feeding  (4)  with  an 
ensemble  of  f/(0,l)  randomly  generated  numbers. 
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Equivalently,  an  ordered  series  is  input  to  (8)  having  the  form 


Construction  Series  = 


N-l 
N  1 


(9) 


Previous  research  established  that  using  a  readily  available  uniform  f/(0,l) 
generator  will  not  produce  a  stable  oc—P  set  for  a  given  /  [ll].  Use  of  a  uniform  random 
number  generator  at  this  stage  introduced  unacceptable  variability  to  the  process,  making 
it  difficult  to  determine  the  optimal  a-p-y  set.  Equation  (9)  can  be  thought  of  as  an 
ordered  collection  of  random  values  from  a  "stable"  uniform  C/(0,1)  random  number 
generator.  In  [1 1]  it  was  established  that  a  value  of  N  «  100,000  produces  a  stable  a-p-y 
set.  In  order  to  speed  up  the  computational  process,  N  =  1000  was  successfully  used. 
This  procedural  convenience  did  not  affect  the  desired  characteristics  of  the  o°  generator. 
Testing  of  the  a°  generator  used  in  combination  with  a  uniform  random  number  generator 
converged  to  the  published  mean  a°  and  s(e)  values  when  enough  samples  were 
produced. 

The  resulting  sequence  is  converted  to  decibels  by  using: 

<7°  (dB)  =  20  logio  (<7i)  (10) 

The  sequence  of  values  in  (10)  is  then  averaged  to  yield  the  mean  cf  value  [3] 


ro 


n  i=i 


HdB) 


(11) 
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The  standard  deviation  of  the  ensemble,  s(0),  is  given  by: 

2 

(12) 


s(e)  = 


/=1 _ 


n-\ 


Each  mean  (t°  and  s(0)  value  is  built  into  the  two-dimensional  arrays  according  to 
their  a-fi  values  (Figures  4  and  5).  The  abscissa  and  ordinate  correspond  to  P  and  a, 
respectively.  For  example,  suppose  a  set  of  a-/?-/ values  is  specified:  aj-pj-/]  (keep 
in  mind  that  each  surface  is  built  for  a  fixed  x  value).  Each  resulting  mean  a°  and  s(0) 
values  have  a]-Pj  coordinates  within  their  respective  arrays.  These  surfaces  constitute 
"search  grids"  or  surfaces  where  an  exhaustive  search  seeks  to  determine  an  «-/?pair  that 
yields  the  published  mean  a°  and  s(0)  values  with  minimum  error. 

An  adequate  measure  of  the  error  for  this  process  is  given  by  [1 1] 


A-B 
error  =  — ^ — 


(13) 
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The  value  A  represents  a  proposed  value  for  a  variable  (i.e.,  A  =  or  A  =  s(0))  while  B 
is  the  desired  (or  measured)  value.  A  measure  of  error  is  obtained  for  both,  mean  a”  and 
s(0)  surfaces,  by  subtracting  and  normalizing  these  by  the  mean  (T°  or  s(0)  value  of 
interest: 


tv 


error  -  ABS 


1  N 

i  i 

M  MN 


-B 


B 


(14) 


V  J 

where  the  array  of  dimension  MN  can  be  either  the  a°  or  s(0)  array  (surface).  B  is  the 
desired  a°  or  s(0)  value,  while  ABS  stands  for  absolute  value. 

The  total  error  (TE)  is  given  by 

TE  =  cf  error  +  s(  0)  error  (15) 

where  (f  error  and  s(G)  error  are  calculated  using  Equation  (14).  The  TE  array  is  a 
surface  of  the  error  incurred  given  the  calculated  and  desired  mean  a°  and  s(0)  values 
(Figure  6).  The  search  algorithm  finds  the  minimum  combined  error  and  determines 
which  a-fi  pair  produced  this  minimum  value.  Notice  the  minimum  error  combination  at 
P  ~  I  and  a  ~  0.5  (Figure  6).  The  algorithm  further  refines  the  estimate  until  the  desired 
level  of  decimal  accuracy  is  reached.  The  algorithm  will  also  stop  after  it  exceeds  a  pre¬ 
set  threshold  and  a  satisfactory  a-P  pair  has  not  been  found. 
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The  criterion  used  to  select  the  range 
of  a-p-y  values  is  based  on  experimental 
experience.  In  this  research,  it  was  found 
that  by  extending  the  a  sequence  to  include 
values  smaller  than  0.16  (i.e.,  0.06) 

convergence  was  possible  for  all 
distributions. 

In  [11],  it  was  found  that  a  search  grid  that  is  dense  for  small  values  of  a-p  and 
sparse  for  large  values  of  these  variables  provides  an  adequate  search  surface. 
Alternatively,  a  "partial  surface"  approach  was  tested.  In  this  method,  only  those  areas  in 
the  surfaces  that  contained  the  mean  and  s(0)  values  of  interest  were  kept.  However, 
this  method  failed  to  converge  to  a-p-y  sets  that  produced  minimum  error  combinations 
for  the  values  of  interest.  The  variable  density  search  method  was  used  instead  with  the 
expected  results.  The  use  of  variable  density  search  spaces  results  in  computational 
timesaving. 

Properties  of  g”  Generating  Function 

Equation  (8)  has  a  discrete  output  dependent  on  the  combinations  and  the 

ensemble  of  C/(0,1)  random  values.  A  histogram  of  the  output  indicates  that  the  resulting 
distributions  can  be  regarded  as  Weibull  (Figure  7).  A  percentile  plot  suggests  this  to  be 
true  (Figure  8),  which  is  representative  of  most  generated  distributions  [3].  The  a-p-y 
values  determine  the  shape  of  the  function's  histogram.  The  value  of  y  can  be  increased 
or  decreased  to  adjust  the  shape  of  the  generated  histogram  and  conform  it  to  the 


Figure  6.  Total  Error  Array 
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Output  Histogram  for  Shmbs.  S*band,  HH  polarization.  6=  8D 


-1.9 


Percentile  Plot  for  Shrubs.  S-band,  HH  polarization,  6  =  80 


Figure  7.  a°  Generator  Output  Histogram.  Figure  8.  Percentile  Plot  of  <7°  Generator 

Output  in  Figure  7. 


measured  data  histogram.  Adjustments  to  y  have  little  effect  when  the  generated 
histogram's  shape  departs  significantly  from  the  measured  data  histogram.  The  bin  width 
of  the  histograms  is  0.5  dB. 

Data  Quality  and  Pre-Processing 

The  data  used  in  this  research  represents  a  calibrated  compilation  of  a° 
measurements  from  several  studies  [8:79].  It  is  the  most  comprehensive  to  date,  although 
data  quantity  varies  by  terrain  type  and  measurement  conditions. 

The  data  sets  were  qualified' by  evaluating  whether  the  mean  o°  and  s(0)  data 
points  could  be  interpolated  to  yield  "reasonably"  shaped  curves  (such  as  seen  in  Figure 
3).  A  data  set  is  any  table  entry  such  as  "1-5,  X  Band,  HH  Polarization"  from  the 
reference  database  [8:91].  The  number  of  source  data  points  (N)  supporting  each  mean 
<5°  value  and  the  shape  of  the  interpolated  curves  were  also  considered.  This  data  pre¬ 
processing  was  necessary  due  to  the  high  variability  in  N  within  and  across  data  sets 
(Figure  9)  [8]. 
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Figure  9.  Partition  of  Data  Points  by  Even  Figure  10.  Partition  of  Data  Points  as  in 
and  Odd  Incidence  Angles.  Previous  Figure. 


The  mean  c°  curve's  shape  is  strongly  linked  to  N  (Figures  9,  10,  and  11).  Note  that  the 
s(6)  curves  also  reflect  the  same  anomalies  as  the  a°  curves  (Figure  12). 

Figure  10  shows  cases  where  the  a"  data  sets  have  markedly  different  N  for  the 
first  six  data  points,  yielding  a  highly  irregular  mean  a°  curve.  This  curve  was  corrected 
by  discarding  the  five  degree  increments  and  raising  the  remaining  g°  data  points  by  an 
amount  such  that  the  resulting  curve  had  the  characteristic  mean  contour.  A  similar 
process  was  applied  to  other  mean  a°  curves  exhibiting  the  same  behavior. 

A  case  where  no  data  point  was  discarded  is  illustrated  in  Figure  11.  Inclusion  of 
the  data  point  6  =  30°  did  not  cause  significant  degradation  of  the  interpolated  mean  a° 
curve. 

Another  pre-processing  example  is  provided  by  Figure  12.  At  first  glance,  the 
data  points  at  0  =  50°  and  55°  seem  to  be  outliers.  However,  these  mean  a°  values  were 
generated  with  a  higher  N  than  the  rest  of  the  points.  One  could  rightly  argue  that  the 
higher  the  number  of  sample  data  points,  the  closer  the  measurement  would  be  to  its  real 
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Band,  HV  Polarization 


value.  However,  it  was  decided  to  discard  these  data  points  in  favor  of  the  remaining 
others.  Keeping  the  majority  of  the  data  points  provided  us  with  more  instances  to  test 
the  generator's  simulation  capabilities.  Moreover,  this  assures  that  the  a°  generator 
output  follows  closely  the  accepted  theory. 

Simulation  of  the  data  collected  over  urban  areas  presented  a  special  situation. 
Analysis  of  the  data  indicated  badly  behaved  mean  g°  curves.  This  is  no  surprise  since 
the  database  authors  warn  readers  that  the  mean  a°  curves  of  urban  areas  may  not  portray 
accurately  the  true  backscattering  response  of  this  terrain  type.  It  is  reasonable  to  infer 
that  these  urban  areas  may  contain  strong  comer  reflectors  (buildings,  wall  layouts,  etc.) 
which  may  influence  measurements  at  certain  radar  frequencies  and  incidence  angles. 
Incidentally,  this  terrain  type  has  the  one  of  the  lowest  number  of  supporting  studies  in 
the  database,  i.e.,  lowest  N.  This  terrain  type  is  included  in  this  research  for  experimental 
purposes  [8:87]. 
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In  general,  the  ten-degree  increments 
were  the  best  documented  in  the  database,  as 
well  as  the  group  to  approximate  closely  the 
theoretical  curve  behavior.  Hence,  these 
were  the  c°  data  source  points  most  often 
used. 

Figure  13.  Spline  Interpolation  ofa°  Data 

Points  for  Short  Vegetation.  Interpolation  of  the  mean  a°  and  s(0) 

curves  was  necessary  to  achieve  the  required  1°  angular  resolution.  The  primary 

reference  applies  an  iterative  least  squares  and  non-linear  curve-fitting  routine  to 

interpolate  these  curves  for  all  incidence  angles  [8].  However,  it  is  desired  to  accurately 

reproduce  the  mean  and  s(0)  measurements  at  each  angle  where  data  is  available. 

Hence,  an  alternative  method  of  interpolation  by  cubic  splines  was  used.  A  spline  curve 

fitting  routine  connects  every  known  data  point  with  a  smooth  curve  that  yields  a  good 

approximation  to  the  theoretical  a°  curves  shown  in  Figure  3.  For  these  purposes, 

MATLAB®'s  spline  routine  was  used.  In  Figure  13  a  "best  case"  interpolation  example  is 

shown  for  short  vegetation  at  L-band,  which  closely  resemble  the  curves  in  Figure  3. 

A  similar  problem  as  the  one  described  above  for  data  points  from  angle  to 

angle  within  a  data  set  occurs  across  data  sets.  In  other  words,  for  the  same  terrain  type, 

radar  band,  and  incidence  angle,  N  generally  differs  for  HH,  HV,  and  VV  polarizations. 

Such  an  instance  is  illustrated  by  shrubs  at  0  =  5”  (L  band),  N  =  204,  74,  and  160  for 

horizontal,  cross,  and  vertical  polarizations,  respectively  [8:200].  This  peculiarity  reflects 

on  the  G°  measurement  history  rather  than  on  the  database's  completeness.  This  may 

explain  why  the  resulting  shape  of  the  interpolated  curves  can  be  irregular. 


Co-Pol  and  Cro-Pol  Components  for  Short  Vegetation/L-Band 
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After  applying  the  analysis  outlined  above,  129  data  sets  out  of  155  were  deemed 
qualified  for  interpolation  and  subsequent  processing. 

Validation  of  Results 

A  multinomial  test  of  proportions  was  implemented  to  test  how  well  the  generated 
histograms  (output  of  Equation  8)  matched  the  measured  (observed)  data  histograms. 
This  is  a  chi-square  goodness-of-fit  test  variation  where  two  distributions  having  a 
different  total  number  of  counts  are  compared.  For  instance,  the  output  of  the  generator 
was  generally  about  100,000  counts.  The  observed  counts  (same  as  N)  varied  from  the 
low  tens  to  approximately  5,000  observations  per  histogram.  To  assess  distribution 
likelihood  between  generated  and  observed  count  histograms,  the  predicted  counts 
histogram  is  normalized.  Then,  each  of  its  bins  is  multiplied  by  the  total  number  of 
observed  counts.  Consequently,  the  sum  of  all  generated  (or  predicted)  counts  equals  the 
total  number  of  observed  counts.  Predicted  and  observed  counts  are  grouped  in 
histograms  having  k  bins  or  cells.  Each  bin  has  width  of  0.5  dB,  consistent  with  the 
database's  convention  [8]. 

Next,  a  null  hypothesis  and  an  alternative  hypothesis  are  proposed.  The  null 
hypothesis  is  that  the  predicted  counts  histogram  is  a  good  fit  to  the  observed  counts 
histogram.  The  alternative  hypothesis  logically  follows  [3:586]. 

The  following  test  statistic  is  used  to  test  the  null  hypothesis: 
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The  test  statistic  value  is  compared  to  the  chi-square  value  obtained  from  tables  or  from  a 
computerized  algorithm.  The  null  hypothesis  is  rejected  when  the  test  statistic  is  larger 
than  the  chi-square  value. 

The  chi-square  value  depends  on  two  factors:  the  level  of  significance  and  the 
number  of  degrees  of  freedom.  The  level  of  significance  can  be  regarded  as  the  tolerance 
of  the  test.  A  level  of  significance  of  0.05  is  used  in  this  research,  which  is  considered  a 
high  level.  In  other  words,  if  the  null  hypothesis  can  not  be  rejected,  we  can  be  95% 
confident  that  this  result  was  not  due  to  chance.  The  number  of  degrees  of  freedom  is 
defined  as  k-l  where  k  is  the  number  of  bins  per  histogram  [3:586]. 

In  a  chi-square  test,  each  bin  must  have  a  minimum  number  of  counts.  In  this 
research  a  minimum  number  of  one  is  used  [5:711].  Due  to  the  normalization  of  the 
generated  data,  the  number  of  predicted  counts  in  any  one  bin  may  be  less  than  one. 
Adjacent  bins  are  combined  together  to  achieve  the  minimum  number  of  counts.  The 
observed  count  bins  are  collapsed  accordingly. 


26 


IV.  Experimental  Results 


This  section  discusses  the  results  and  issues  encountered  in  the  generator 
development.  The  main  focus  in  this  investigation  is  to  develop  the  c°  generator  for  all 
types  of  terrain  and  measurement  conditions  published  in  the  database  [8]. 

Experimental  Data  Validation 


Several  issues  ensued  in  the 
validation  phase.  Perhaps,  the  most 
significant  problem  was  the  mixture  of  well 
and  poorly  documented  incidence  angles 
within  a  data  set.  Although  a  data  set  may 
have  been  qualified  after  pre-processing, 
this  does  not  imply  that  all  of  its  data  points 
were  supported  by  a  high  N  (see 
Chapter  EH,  Data  Quality  and  [8:210,215]). 
Validation  results  were  affected  by  the  value  of  N.  As  N  increased  (hundreds  to 
thousands),  any  small  deviation  between  the  predicted  and  observed  counts  histograms 
resulted  in  rejection  of  valid  null  hypotheses.  This  is  shown  in  Figure  14,  where  TS  is  the 
test  statistic,  NTS  is  the  normalized  test  statistic  (to  be  discussed  shortly),  and  is 

the  chi-square  value. 

On  the  other  hand,  when  the  measured  counts  were  low  (N  <  100)  the  test  did  not 
reject  ambiguous  or  invalid  null  hypotheses  in  almost  all  instances  (Figure  15).  Visual 
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Figure  15.  Grass/S-Band/HV  pol./0  =  5°,  Figure  16.  Chi-square  Test  Success  Rate 
N  =  27,  TS  =  25.  NTS  =  92.  =  3 1 .  Distribution 

inspection  of  the  histograms  indicated  that  a  reasonable  conclusion  could  not  be  reached 

due  to  the  low  number  of  observed  samples  and  sparse  histograms. 

The  chi-square  test  seemed  to  be  prone  to  reject  the  null  hypothesis  when  N  was 

high  and  viceversa.  To  test  this  suspicion,  a  ratio  test  was  devised: 


ratio  test  = 


(17) 


where  the  numerator  is  the  test  statistic  in  Equation  (16)  and  the  denominator  is  the  chi 
square  value  for  a  significance  level  a  and  ^-1  degrees  of  freedom.  When  the  null 
hypothesis  is  not  rejected,  the  numerator  will  be  smaller  than  the  denominator  yielding  a 
ratio  less  than  one.  The  opposite  occurs  when  the  null  hypothesis  is  rejected.  The 
histogram  in  Figure  16  shows  what  fraction  of  chi-square  tests  within  each  bin  yielded  a 
ratio  test  less  than  one.  It  is  seen  that  the  validation  test  tends  not  to  reject  the  null 
hypothesis  when  N  ranges  in  the  low  tens  and  hundreds.  The  test  behavior  was 
incompatible  with  the  goal  of  validation  process  automation.  An  alternative  approach,  a 
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modified  chi-square  test,  was  tested  and  implemented.  This  test  yielded  good  results  for  a 
greater  range  of  N  values  as  confirmed  by  a  detailed  review  of  all  validation  test  results. 

The  test  is  implemented  as  an 
ordinary  chi-square  test  except  that  both, 
predicted  and  observed  bin  counts,  are 
normalized  to  one  hundred.  Due  to  this 
normalization,  the  test  statistic  in  Equation 
16  is  called  the  normalized  test  statistic 
(NTS).  This  test  performed  better  than  the 
regular  chi-square  test,  especially  for  low  N. 
Figure  17  shows  the  ratio  test  results  for  the  normalized  chi-square  test  (Equation  17).  It 
is  seen  that  the  test's  success  rate  is  better  distributed  along  the  range  of  samples  shown. 

However,  the  modified  chi-square  test  did  not  reject  inconsistent  null  hypothesis 
for  certain  high  N  situations  where  the  measured  data  histograms  had  a  wide  distribution 
(Figure  18).  These  were  data  sets  that  the  functional  form  (Equation  8)  could  not 
accurately  represent,  most  frequently  for  low  incidence  angles  in  the  region  0  <  20°. 
Fortunately,  failure  of  the  test  was  extremely  rare  (less  than  5  occasions  in  958). 

On  the  other  hand,  the  modified  chi-square  test  did  fail  generated  distributions 
that,  to  the  naked  eye,  appeared  acceptable  (Figure  18(b)).  This  is  noticeable  especially 
for  distributions  with  100  <  N  <  200,  approximately.  Of  course,  the  test  is  designed  to 
take  into  account  bin  height  differences  between  generated  and  measured  histograms  the 
naked  eye  cannot  detect.  A  list  of  these  distributions  is  included  in  the  appendix. 


Figure  17.  Normalized  Chi-square 
Non-rejection  Rate. 
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Figure  18(a).  Validation  Results  for  Short 
Vegetation:  TS  =  87,  NTS  =  63,  ^  =  65, 

N=535 

The  observations  outlined  above  serve  as  an  introduction  to  a  more  detailed 
treatment  of  the  results  by  terrain  type.  Figures  19-27  show  graphical  output  of  the 
generator  (dark  line)  on  top  of  measured  data  histograms  (rectangular  bins).  The  title  of 
each  figure  describes  the  terrain  type,  radar  band,  polarization,  and  angle  of  incidence. 
For  example,  aLHHO  means  soils  and  rocks,  L  band,  HH  polarization,  and  0  =  0°.  The 
test  statistic,  chi-square  value,  and  normalized  test  statistic  values  are  in  the  appendix. 
Soils  and  Rocks 

In  general,  this  terrain  type  (type  "a")  is  well  documented  with  N's  ranging  in  the 
thousands  for  some  cases  [8:130].  The  main  discrepancies  between  the  generator  and 
the  measured  data  were  found  mostly  at  low  incidence  angles  (0°  <  0  <  30°).  Excellent 
agreement  was  observed  for  L,  S,  and  Ku  band  for  N  >  100.  C  band  performance  was 
acceptable  (close  to  70%  success  rate)  while  X  band  performance  was  poor 
(Figures  19(a-i)). 


Figure  18(b).  Validation  Results  for 
Grasses:  TS  =115,  NTS  =  67,  ^  =  34, 
N=  172 


30 


1)  distribution  rejection  for  low  incidence  angles,  2)  rejection  due  to  low  number  of 
distribution  samples  of  fair  distributions,  and  3)  rejection  of  fair  distributions.  About 
97%  of  all  distributions  have  N  <  100.  The  generator  performance  was  poor  for  this 
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terrain  type  (Figures  20(a-i)).  Special  distribution  instances,  such  as  bimodal 
distributions,  were  noted  in  this  terrain  type.  Figure  20  (h)  illustrates  a  bimodal 
distribution  with  N  =  70.  It  is  reasonable  to  suspect  that  this  situation  is  due  to  low  N. 


Figure  20  (a)  bXHHO 
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Figure  20  (d)  bXHVSO 
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Figure  20  (g)  bKuHVbO 


«• . 1 . :• . r 

ge- . J . ; . i 


<r  m 

Figure  20  (b)  bXHH60 
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Figure  20  (h)  bKuHH40 
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Figure  20  (c)  bXHVO 
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Figure  20  (f)  bKuHVO 
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Grasses 


This  terrain  type  (type  "cl")  featured  high  N  for  many  data  sets.  The  generator 
performed  well  for  all  radar  bands  except  L  band  (Figure  21  (a-i)). 
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Figure  21  (g)  clKuHVSO 
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Figure  21  (b)  clLHHSO 
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Figure  21  (e)  clCHH60 
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Figure  21  (h)  clKuVVbO 
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Shrubs 


The  a°  generator  performed  satisfactorily  for  this  terrain  type  (type  "c3")  for  most 
incidence  angles  (Figures  22  (a-i)).  Discrepancies  were  noted  in  the  region  0°  <  0  <  30® 
for  all  data  sets,  especially  when  100  <  N  <  200.  The  generator  had  excellent  to 
acceptable  performance  for  all  bands.  The  average  success  rate  was  65%. 
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Figure  22  (d)  c3SHV60 
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Figure  22  (f)  c3XHV40 
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Short  Vegetation 


This  terrain  type  (type  "c")  featured  good  agreement  for  most  well  documented 
data  sets  (N  >  100)  for  0  =  30°  and  above  (Figures  23  (a-g)).  Average  success  rate  was 
about  80%.  Excellent  to  good  performance  for  all  radar  bands.  Incidentally,  this  terrain 
type  features  very  large  N  [8:247]. 
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Figure  23  (d)  cSHH80 
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Figure  23  (f)  cKaHVSO 
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Roads 


This  terrain  type  (type  "d")  was  characterized  by  N  <  100  and  sparse  histograms. 
Some  of  the  distributions  were  well  suited  for  description  using  uniform  distributions 
(Figures  24  (a-c)).  Consequently,  the  generator  performance  was  poor.  However,  a  fair 
assessment  of  the  modeling  capabilities  of  the  a°  generator  for  this  terrain  type  requires  a 
higher  N. 


Figure  24  (a)  dXHH20  Figure  24  (b)  dKuHH60  Figure  25  (c)  dKaVV  10 

Dry  Snow 

For  this  terrain  type  (type  "f '),  the  a"  generator  performance  varied  from  excellent 
to  good  for  distributions  with  N  >  100  and  L,  S,  C,  and  Ku  radar  bands  (Figures  25  (a-i)). 
The  success  rate  was  about  70%  for  these  bands.  Also,  visual  inspection  of  rejected 
histograms  demonstrated  good  agreement  between  the  generator  and  measured  data 
(appendix). 


Figure  25  (a)  fLHV20 
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Figure  25  (g)  fKuVVTO  Figure  25  (h)  fKaHH50  Figure  25  (i)  fWVVlO 


Wet  Snow 

For  wet  snow  (type  "g"),  the  generator  functional  form  yielded  good  agreement 
for  C,  X,  Ku,  and  Ka  bands  with  a  success  rate  of  about  75%  (Figures  26  (a-i)). 


Performance  was  poor  for  L  and  S  bands  at  low  incidence  angles. 


Figure  26  (a)  gLHH70  Figure  26  (b)  gLHV70  Figure  26  (c)  gSHHSO 
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Figure  26  (d)  gSHVSO 
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Figure  26  (g)  gXHVSO 


9SVV50 


Figure  26  (e)  gSVVSO 

gKlVVTD 


Figure  26  (h)  gKuVVTO 
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Summary  of  Results 


The  experimental  results  demonstrate  the  soundness  of  the  <7**  generator  approach. 
The  function  in  (8)  has  broad  applicability  to  several  terrain  types,  wave  polarization,  and 
incidence  angle  accounting  for  an  overall  success  rate  of  72%  for  distributions  with 
source  data  points  (N)  greater  than  or  equal  to  100.  These  observations  are  documented 
in  the  appendix. 

In  general,  the  o°  generator  performed  well  for  short  vegetation,  shrubs,  and  soils 
and  rocks  except  in  the  sector  0°  <  0  <  20°.  The  distributions  at  these  incidence  angles 
exhibited  strong  Rayleigh  behavior.  For  wet  snow  and  dry  snow,  good  agreement  could 
be  observed  for  all  incidence  angles.  Roads  and  trees  have  the  lowest  N  of  all  t5q)es; 
hence,  rejection  rates  were  high. 

It  is  difficult  to  explain  why  the  success  rate  is  not  higher  based  on  the  limited 
amount  of  data  and  its  distribution  across  terrain  types.  For  instance,  the  a°  generator 
performed  very  well  for  short  vegetation  and  shrubs  except  for  low  incidence  angles. 
Incidentally,  these  terrain  types  have  the  highest  total  N.  The  opposite  occurs  for  terrain 
types  such  as  trees  and  roads  (appendix).  It  would  be  speculative  to  assert  that  as  N  is 
increased  for  all  terrain  types,  the  generator  will  perform  at  least  as  good  as  for  shrubs 
and  short  vegetation.  However,  it  is  a  possibility  given  the  experimental  results. 

The  modified  chi-square  test  proved  to  be  a  proper  metric.  However,  when  the 
normalized  test  statistic  and  the  chi-square  value  are  commensurate,  it  is  advisable  to 
visually  judge  the  generated  and  measured  histogram  likelihood.  This  is  also  true  when  N 
<  200,  approximately,  and  when  observed  data  histograms  have  wide  distributions. 
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V.  Conclusions  and  Recommendations 

The  functional  form  in  (8)  accurately  converged  to  the  published  mean  g°  and  s(9) 
values  for  nine  terrain  types,  radar  bands,  the  linear  pair  of  wave  polarizations,  and 
incidence  angles  lower  than  80°.  Equation  (8)  was  also  shown  to  exhibit  the  same 
frequency  of  occurrence  histograms  for  approximately  72%  of  measured  a°  distributions 
with  database-source-points  (N)  >  100.  Above  average  results  were  noted  for  short 
vegetation  (78%),  soils  and  rocks  (75%),  and  grasses  (74%).  However,  experimental 
results  suggest  that  an  alternate  functional  form  (or  modification  to  the  current  form)  may 
be  necessary  to  accurately  portray  short  vegetation,  soils  and  rocks,  grasses,  and  shrubs 
for  low  incidence  angles.  This  inconsistency  was  not  observed  for  dry  and  wet  snow 
where  the  ct°  generator  performed  well  for  all  incidence  angles  (72%  and  69%  success 
rates,  respectively).  Roads  and  trees  require  a  higher  N  in  order  to  establish  the  c° 
generator  suitability.  Terrain  type  “urban”  was  not  validated  since  all  of  its  distributions 
exhibited  low  N. 

A  consequence  of  the  variability  in  N  is  that  within  c°  data  sets  that  had  a  high 
number  of  failed  distribution  tests,  a  success  could  be  found.  The  converse  was  also  true. 
This  variation  affected  the  reliability  of  the  validation  tests.  Certainly,  data  quantity  and 
quality  is  an  important  factor  to  be  considered  for  future  work. 

The  following  recommendations  are  given: 

(1)  Campaigns  to  measure  ct°  for  several  terrain  types  to  avoid  the  variability  in  N 
found  within  data  sets  for  all  incidence  angles.  This  will  help  to  precisely  evaluate  the  a° 
generator  in  a  “controlled”  setting. 
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(2)  Explore  alternative  functional  forms  that  resemble  the  behavior  of  a°  curves. 


Available  statistical  software  packages  can  be  used  for  this  purpose. 

(3)  A  modification  to  the  functional  form  to  appropriately  model  the  (f  from 
certain  terrain  types  for  low  incidence  angles. 

(4)  Further  refinement  of  the  normalized  chi-square  test  to  signal  when  the 
computational  process  achieves  the  desired  g°,  s(6)  values,  and  frequency  of  occurrence 
histogram  simultaneously. 
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Appendix:  Simulated  Data  Sets 


This  appendix  contains  the  generated  data  for  this  research.  A  summary  for  all 
generated  and  validated  distributions  is  included  in  page  44.  The  data  set  groups  included 
are  udmodO,  udmodl,  udmodl,  udmod3,  and  udmod4a. 

The  name  of  the  data  set  group  denotes  the  data  set  composition.  For  example, 
udmodO  includes  all  published  data  sets  whose  mean  values  were  supported  by  N  >  16. 
Hence,  zero  distributions  had  N  <  16.  The  data  set  group  udmodda  is  composed  of  data 
sets  that  have  four  or  more  incidence  angles  with  N  <  16. 

Each  data  set  table  in  the  appendix  has  14  labeled  columns; 


(1)  lA  =  incidence  angle 

(8)  Gamma  =  gama  parameter 

(2)  Mean  =  published  mean  <f  value 

(9)  TstStat  =  test  statistic 

(3)  ResMean  =  generated  mean  a°  value 

(10)  OAChi  =  chi-square  value 

(4)  StD  =  published  s(0) 

(11)  NormTS  =  normalized  test 

statistic 

(5)  ResStD  =  generated  s(0) 

(12)  N  =  source  data  points 

(6)  Alpha  =  alpha  parameter 

(13)  NTS/OACHI  =  ratio  test  for 

NormTS 

(7)  Beta  =  beta  parameter 

(14)  TS/OACHI  =  ratio  test  for 

TstStat 
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A  value  of  -1  in  column  (10)  indicates  that  the  number  of  histogram  bins  in  a 
distribution  was  less  than  two.  Validation  tests  were  not  done  on  distributions  with  less 
than  two  bins. 

The  following  line  in  the  table  describes  the  terrain  type  and  measurement 
conditions.  For  instance,  ACHHO  denotes  terrain  type  A  (also  referred  as  a),  C-Band,  HH 
polarization,  and  y  value.  The  terrain  type  code  in  italics  is  also  used  in  this  research, 
specifically  when  validating  the  generated  distributions.  The  following  list  describes  all 
terrain  types,  radar  bands,  and  polarizations  [8:88]: 

(1)  A,  a  =  soils  and  rocks  (10)  L  =  1-2  GHz 

(2)  B,  Zj  =  trees  (1 1)  S  =  2-4  GHz 

(3)  C,  cl  =  grasses  (12)  C  =  4-8  GHz 

(4)  D,  c5  =  shrubs  (13)  X  =  8-12  GHz 

(5)  E,  c  =  short  vegetation  (14)  Ku  =  12-18  GHz 

(6)  F,  (i  =  roads  (15)  Ka  =  30-40  GHz 

(7)  G,  e  =  urban  (16)  W  =  90-100  GHz 

(8)  H,/=  dry  snow  (17)  HH,  VV  =  co-polarized  components 

(9)  I,  g  =  wet  snow  (18)  HV  =  cross-polarized  components 

The  summary  table  gives  a  bird's  eye  view  of  validation  results  by  terrain  type  and 
radar  band.  Validation  tests  are  for  normalized  and  regular  chi  square  tests  for  two  range 
of  samples:  1)  N  >  16  and  2)  N  >  100.  The  label  GE  and  LE  mean  equal  or  greater  than 
and  equal  or  less  than,  respectively. 
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Added  Histograms  (N  >=  100) 

aCHH30/20/I0 

aCHV50/10 

C1CVV40  (picture)  TS=67  CS=34 

clKuHHSO 

clKuVV20 

clSHHO 

clSHVO  (picture)  TS=54  CS=47 

clSVV30,40 

C3CHV20 

C3CVV40/50 

C3XHH40 

C3XHV20 

C3XVV30/40 

C3SHH40 

C3SVV30/40 

c3KuHH30 

c3KuVV20 

C3XHV20 

C3XHH40 

CSVV40  (picture)  TS=52  CS=40 

CSHV60 

CCVV40 

cKuVV20 

fKuHHSO 

fKuVVSO 

fKaHHSO 

fLHVO  (picture) 

fSHV50 

fSVV50,30,20 

fCVVO 

gCHV70 

gCVV70 

gSHH20,  50  (picture)  TS=43  CS=41 

gCHV30 

gLHHSO 

gLVVSO 

Subtracted 

clCVV  10  (picture)  TS=77  CS=52 
clKuHHlO 

cCHHO  (picture)  TS=57  CS=67 
cSHHO 


Total  Samples  per  Tej  rjin  Type 

soils  34291  KM 

trees  3269  I'i 

grass  43013  12  4Mo 

shrubs  107959  3  % 
short  veg.  134648  2  9% 

roads  1195  O.J'/o 

dry  snow  11880  3  3% 

wet  snow  10341  :% 

Some  "success"  rates  :  or  N  >=  100 

short  veg.  78% 

soils  75% 

grasses  74  % 

dry  snow  72% 

'wet  snow  69% 
shrubs  66  % 
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udmodO 


lA  Mean  ResMean  StD  ResStD  Alpha  Beta  Gamma  TstStat  OAChi  NormTS  N  NTS/  TS/ 
at  hho  OACHI  OACHI 
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lA 

AL 
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WO 
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lA 

AS 

Mean 

HH2 
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ll^KE 
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IB 
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3.93 

0.63 
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udmodO 


lA  Mean  ResMean  StD  ResStD  Alpha  Beta  Gamma  TstStat  OAChi  NormTS  N  NTS/  TS/ 

AS  WO  _ OACHI  OACHI 


m 

■s 

7.6 

0.68 

1.95 

0 
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73.3 

153.6 
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2.10 

5.53 

n 

-0.9 

-0.9 

5.8 
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43.1 
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1.33 

■Q 
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mm 

4.1 

0.426 
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0 
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60.5 

19.9 
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0.33 

1.45 

■B 

Kl 

mm 

4.7 

0.426 

0.611 

0 

72.4 

60.5 

20 

362 

0.33 

1.20 

-9 

-9 

13 

4.6 

0.417 

0.451 

0 
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61.7 

29.2 

442 

0.47 

2.09 

Rn 

-12.6 

-12.6 

5.2 

5.2 
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0.307 

0 
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58.1 

100.9 
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4.57 

m 
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■El 
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0 
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82 
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0.72 

m 

-17.8 
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3 
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0 
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lA  Mean  ResMean  StD  ResStD  Alpha  Beta  Gamma  TstStat  OAChi  NormTS  N  NTS/  TS/ 

BK  UHH5  _ OACHI  OACHI 
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0.126 
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lA  Mean  ResMean  StD  ResStD  Alpha  Beta  Gamma  TstStat  OAChi  NormTS  N  NTS/  TS/ 

BK  UHV5  _ OACHI  OACHI 
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0.05 

45.2 

38.9 

70.7 

64 

1.82 

1.16 

-12.2 

-12.2 

3.3 

3.3 

0.408 
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lA  Mean  ResMean  StD  ResStD  Alpha  Beta  Gamma  TstStat  OAChi  NormTS  N  NTS/  TS/ 

BK  UVV5  _ OACHI  OACHI 
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lA  Mean  ResMean  StD  ResStD  Alpha  Beta  Gamma  TstStat  OAChi  NormTS  N  NTS/  TS/ 

BX  HV5  _ OACHI  OACHI 
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